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Abstract — This work focuses on a large-scale multi-cell multi¬ 
user MIMO system in which L hase stations (BSs) of N antennas 
each communicate with K single-antenna user equipments. We 
consider the design of the linear precoder that minimizes the 
total power consumption while ensuring target user rates. Three 
configurations with different degrees of cooperation among BSs 
are considered: the coordinated beamforming scheme (only 
channel state information is shared among BSs), the coordinated 
multipoint MIMO processing technology or network MIMO 
(channel state and data cooperation), and a single cell beam¬ 
forming scheme (only local channel state Information is used 
for beamforming while channel state cooperation is needed for 
power allocation). The analysis is conducted assuming that N 
and K grow large with a non trivial ratio K/N and imperfect 
channel state information (modeled by the generic Gauss-Markov 
formulation form) is available at the BSs. Tools of random 
matrix theory are used to compute, in explicit form, deterministic 
approximations for: (i) the parameters of the optimal precoder; 
(ii) the powers needed to ensure target rates; and (Hi) the total 
transmit power. These results are instrumental to get further 
insight into the structure of the optimal precoders and also to 
reduce the implementation complexity in large-scale networks. 
Numerical results are used to validate the asymptotic analysis 
in the finite system regime and to make comparisons among the 
different configurations. 


I. Introduction 

The road forward for satisfying the increasing number of 
users and high rate expectations in 5G systems is very high 
spatial utilization. Among the different technologies in this 
context, massive MIMO is considered as one of the most 
promising [l]-[4]. Under the assumption of uncorrelated chan¬ 
nels, if the number of base station (BS) antennas N goes to 
infinity and the number of user equipment terminals (UE) K is 
maintained fixed, the performance of massive MIMO systems 
becomes limited only by the so-called pilot contamination and 
simple matched filter and maximum ratio transmission with 
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no cooperation among cells can entirely eliminate the uplink 
and downlink multicell interference. However, the maximum 
number of antennas at each BS is limited in practice. In such 
a case, interference-aware precoder designs with cooperation 
among cells should be applied for optimal handling of the 
remaining intercell interference. 

Coordinated multi-cell resource allocation is generally for¬ 
mulated as an optimization problem in which a desired net¬ 
work utility is maximized subject to some quality-of-service 
requirements. In this work, we focus on the problem of 
designing the optimal linear precoder in multi-cell networks 
for minimizing the total transmit power while ensuring a set 
of target user rates [5], [6]. This problem has received great 
attention in the last years [7] and is gaining renewed interest 
nowadays due to the emerging research area of green multi 
cellular networks [8]. We specifically consider the downlink 
of a multi-cell multi-user MIMO system in which L BSs 
equipped with N antennas each communicate with K single¬ 
antenna UEs. Within this setting, several configurations with 
different degrees of cooperation can be envisioned [9]. In this 
work, the following three are considered; (/) the coordinated 
beamforming (CoBE) scheme in [10] in which each BS sends 
data to its own users only but channel state information 
(CSI) is shared between the L BSs so that each BS can 
exploit its excess number of spatial dimensions to mitigate the 
interference generated in other cells; this has the advantage 
of not having to distribute all users’ data to all BSs; (ii) 
the fully cooperative scheme, widely known in the literature 
as coordinated multipoint MIMO (CoMP) or also network 
MIMO [9], in which the BSs share the CSI of the UEs 
as well as their data through backhaul links'; and (Hi) a 
single cell beamforming (ScBE) scheme in which each BS 
computes the beamformer only on the basis of the CSI of 
its own UEs while the power allocation requires the BSs to 
share CSI so that the rate constraints can be jointly satisfied. 
In all the above scenarios, under the assumption of perfect 
CSI, the optimal linear precoder is known to be a function 
of some Lagrange multipliers, the computation of which can 
be performed using convex optimization tools or solving a 
fixed-point problem [7]. Although numerically feasible, both 
approaches do not provide any insight into the structure of 
the optimal precoder. Moreover, the computation must be 
performed for any new realization of propagation channels, 
which is too computationally cumbersome when the network 
size becomes large (as envisioned in 5G networks). 

'Although CoMP might refer to a more general setting in the 3GPP 
standard, in this work the term is used to refer to network MIMO. 
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To overcome these issues, we follow the same approach as 
in other works for single- or multi-cell networks [11]-[16] 
and resort to the asymptotic regime where N and K grow 
large with a non trivial ratio K/N. The design and analysis 
of the considered networks is performed under the assumption 
of imperfect CSI (modeled by the generic Gauss-Markov 
formulation form) for the UEs. Unlike most previous works 
[13], [16], the asymptotically optimal values of the Lagrange 
multipliers are computed using recent results from random 
matrix theory [17], which provide us with a much simpler 
means to overcome the technical difficulties arising with the 
application of standard random matrix theory tools. These re¬ 
sults are then exploited to compute explicit expressions for the 
asymptotic signal-to-interference-plus-noise ratios (SINRs), 
which are eventually used to obtain the asymptotic powers 
needed to ensure target rates as well as the asymptotic total 
transmit power. As shall be seen, all the aforementioned 
deterministic approximations are found to depend only on 
the long-term channel attenuations of the UEs, the relative 
strength of interference among BSs, the target rates and the 
quality of the channel estimates. As a notable outcome of 
this work, the above analysis provides a unified framework 
that can be used to compare the considered networks under 
different settings and to eventually get insights on how the 
different parameters affect their performance. Moreover, in the 
same spirit of [18], [19], the provided results can be used to 
derive optimal distributed algorithms that rely only on the ex¬ 
change, among nearby BSs, of long-term fading components. 
Numerical results are used to show that the asymptotic analysis 
well approximate the network performance in the finite system 
regime. 

The main literature related to this work is represented by 

[12] -[16], [20]. Specifically, a single-cell setting is consid¬ 
ered in [12], [20] while a CoBE network is investigated in 

[13] . Unlike [13], we provide closed-form expressions for the 
Lagrange multipliers, which are instrumental to also compute 
closed-form expressions for SINRs and transmit powers. In 
[15], the authors focus on the sum rate of a CoMP under 
the assumption of regularized zero-forcing precoding. In [16], 
the authors provide an asymptotic analysis of all considered 
network configurations but for the simplest case in which 
only two cells are present and CSI is perfect. The analysis 
is also conducted under the restrictive assumption that the 
channels among all UEs within a given cell and the interfering 
BS can be modelled as ^ CA((0, cIat) with e controling 
the interference level between neighbor cells. Moreover, the 
same rate is imposed to all UEs. We importantly show in 
the present article that, within our framework, there is no 
substantial additional difficulty in treating the more general 
setting of interest here which, unlike [16], considers multiple 
cells, different rate constraints, imperfect CSI and models the 
interference between neighbor cells according to the large 
scale fading coefficients. 

The remainder of this paper is organized as follows. Next 
section describes the signal model and revisits the optimal 
linear precoder for the different network configurations. Sec¬ 
tion III deals with the large system analysis. In Section IV, 
numerical results are used to validate the theoretical analysis 


and to make comparisons among the different beamforming 
schemes. Einally, the major conclusions and implications are 
drawn in Section V. All the technical proofs are presented in 
the Appendices. 

Notation: The following notation is used throughout the 
paper. Scalars are denoted by lower case letters whereas bold¬ 
face lower (upper) case letters are used for vectors (matrices). 
The superscripts ^ and ^ denote transpose and conjugate 
transpose. We denote by the identity matrix of order N 
and call On and the dimensional all-zero and all- 
one vectors, respectively. A random vector x ~ CA/^(m, C) 
is complex Gaussian distributed with mean m and covariance 
matrix C. The notation tr(A) stands of the trace of matrix 
A whereas A = diag{a;i,..., denotes a diagonal matrix 
of order N. We use [-ji ^ to denote the {i, fc)th element of 
the enclosed matrix and 0 to indicate the Kronecker product. 
We denote a„ x bn the equivalence relation an — bn ^ 0 as 
n —> oo for two infinite sequences a„ and bn- 

II. System Model 

Consider the downlink of a multi-cell multi-user MIMO 
system composed of L cells, the BS of each cell comprising 
N antennas to communicate with K single-antenna UEs. As 
mentioned previously, we consider three different configura¬ 
tions with different degrees of cooperation: i) the coordinated 
beamforming scheme in [10]; ii) the coordinated multipoint 
processing (or network MIMO) [9]; and in) a single cell 
beamforming scheme. In all these scenarios, we are interested 
in minimizing the total transmit power Pr while satisfying rate 
constraints at the UEs. Under the assumption of perfect CSI at 
the BSs, this problem can be solved using different approaches 
based for example on second order conic programming and 
standard decomposition techniques [7]. Next, we consider a 
finite size system and review the optimal precoder structure 
for the aforementioned schemes. In doing so, we assume that 
the feasibility conditions are satisfied [5]-[7], [13]. 

A. Coordinated Beamforming 

In the CoBE setting, each UE is attached to a specific 
serving BS while receiving interfering data from other BSs. 
As such, we shall use a double index notation to refer to each 
UE e.g., “user k in cell f\ Under this convention, let us thus 
define hijk € as the channel from BS I to UE k in cell j, 
given by hijk = yjdijk'^ijk where is the small- 

scale fading channel assumed to be Gaussian with zero mean 
and unit covariance, i.e., ~ CN{0n,^n), and dijk accounts 
for the corresponding large scale channel fading or path loss 
(from BS I to UE k in cell j). Denoting by gjk € the 
precoding vector intended to UE k in cell j, its received signal 
can be written as 

K 

Vjk = ^jjkSjkSjk + 

L K 

+ Y. iSli+Tljk ( 1 ) 

l=l,l^j i=l 
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with sii G C being the signal intended to user i in cell /, 
assumed independent across {l,i) pairs, of zero mean and 
unit variance, and rijk ~ CAf{0, ct^). Under the assumption of 
target UE rates {rjkiVj, k}, the power minimization problem 
for CoBF can be formulated as: 


min 

{gjfc} 


L K 




( 2 ) 


_ IhgfcgjfcP _ 

1=1,i^3 l = l,l^j2=l 


> Ijk Vj, k 


where 7 ^^ = 2’’ 3 '“ — 1 denote the corresponding SINK 
constraints. Upon existence [7], the unique solution of (2) is 


- 


lil. 


gjfc — \l N ||v7|| Efe ~ VZ^i=lZ^i=l N 

I at) Vj, k where {A*^/7V} are the Lagrange multipliers 

associated to the SINR constraints and are obtained as the 
unique fixed point solution of the following set of equations 
[5]-[7]: 

^-l 

Vj,fc. (3) 


~ V Z^i=l tv “j***Ui) 


X* _ 


(1 + l/ 7 jfc) 


^jjk 


E E 

1=1 i=l 


-1 

hjjk 


The optimal power values {p*/,} such that the SINR 
constraints in (2) are all satisfied with equality. This amounts 
to computing the unique solution of the following set of 
equations [7]: 


1 Pjk 

Ijk N ||v*J|2 


K 

Pji 





L K 


- Y Y 

1=1,l^j i=l 


PH_ 

N 


hyfeV* 


0 Vj, k. 


(4) 


Observe also that the Lagrange multipliers and 

vectors {v*^,} may be thought of as the solution to the 
following dual UL power minimization problem [16]: 


L K 


A 


j=l k=l 


min 

{vjk,Yk} ■“ ■“ ^ 


s.t. 






E + NIn 


(5) 

> Ijk Vj, k 


ifjk 


L(fvv(i.fc) 

with {vjfe} being the receive beamforming vectors. 


B. Coordinated Multipoint Processing 

In the CoMP (or network MIMO) setting, each UE is jointly 
served by all BSs. In other words, there exists no cell-user 
association and thus the UEs can be indexed as k from 1 to 
KL. Let us then denote = [h^^,..., y with hjk G 
the channel from BS j to user k given by h.jk = y^djk'Wjk 
where Wjk G is the small-scale fading and djk accounts 
for the path loss (from BS j to UE k). Denoting by gk G 


the joint precoding vector for UE k, its received signal can be 
written as 

KL 

Vk = hf gfeSfc -f Y hf -f rik (6) 


with Si G C being the signal intended to user i, independent 
across i, of zero mean and unit variance, and Uk ~ CAf{0, a^). 
In the above setting, the power minimization problem takes the 
form: 


min 

{gfc} 


KL 

Y^kSk 


S.t. 


IhfgfcP 

KL 

E |hfg.P+a2 


> Ik 


Vfc 


(7) 


where 7 fc = 2'’'= — 1 with the rate constraint of UE 
k. The solution to (7) is gj = with = 

(Eil^ + Kl) where {A*/(iVL)} are such 

that [5]-[7]: 


Afc = 


(1 + l/7fc)' 


KL 


hf EA*h,hf+iVLI 


NL 


i=l 


V/C. (8) 


As before, the optimal {p^} are computed such that the SINR 
constraints in (7) are satisfied with equality [7]. Then, we 
obtain 


1 pk \hyy ^ Pi |hfv*|^ 

IkNL ||v*f llEf 


0 Vfc. (9) 


As for CoBF, {Xl/{NL)} and {v^} can be obtained as the 
solution to the following dual UL problem [16]: 


KL . 

E Xk 2 

'TTTO' 

NL 

k—l 


S.t. 


Afelvfhfcl 


> Ik 


vf E A,h,hf-k iViljvL vfc 


,i^k 




m 


Vfc. 


C. Single Cell Beamforming 

Inspired by a single cell beamforming scheme (see for 
example [12]), we also consider the following power mini¬ 
mization problem: 


K 


S.t. 


l^fjksjk 


K 


E Ih^feSv 


K 


( 11 ) 

> 7jfc Vfc. 


E Elh^fcgi*P + 


i=i,iGj i=i 

Upon existence, the solution to (11) is given by g*^ = 


with y = + I 


N llv 


jk I 


N 
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where the scalars {A*^/-/V} are such that [5]-[7]: 

(l + l/7ife) 


^jk - 


K 


+iVI 


i-N 


i=l 


yj,k. ( 12 ) 


jjk 


The powers {p*^} are computed such that the SINR con¬ 
straints in (4) are all satisfied [7]. As for CoBF, this requires 
CSI to be exchanged among BSs. This makes the consid¬ 
ered ScBF scheme different from a “classical” single cell 
processing technique in which each BS does not exchange 
any information with the other BSs.^ The reason why it is 
referred to as a single cell scheme is due to the fact that, 
unlike CoBF, the computation of the beamforming vectors 
Vjfc = (E*=i -h lNy\jk of cell j does not 

require knowledge of the fast fading channels {hj^; ^ j, i} 

from all UEs in the other cells to cell j. More details on this 
will be given in Section III-C. 

Similar to CoBF and CoMP, {Xj/^/N} and {v*^,} may be 
thought of as the solution to the following dual UL power 
minimization problem for cell j: 


K 

min y 


Xjk 2 


(13) 


s.t. 






K 




■N 




> ijk Vfc 




with {vjk} being the receive beamforming vectors and = 


III. Large System Analysis 

Let A* and p* denote for the three settings above the 
vectors collecting the Lagrange multipliers and power values, 
respectively. As shown previously, the precoding vectors are 
parameterized by A* and p*, where A* needs to be eval¬ 
uated by solving a set of fixed-point equations. This is a 
computationally demanding task when N and K are large 
since the matrix inversion operation in (3), (8) or (12) must 
be computed several times, with complexity proportional to 
N^KL or {NL)'^KL. Besides, from a practical standpoint, a 
close inspection of the SINR constraints in (4) and (9) reveals 
that the evaluation of p* requires knowledge of all channels 
{hyfc} and {hi}, thus implying some channel exchange proce¬ 
dure within the network at the rate of the fast fading channel 
evolution. This makes the implementation of all the above 
solutions a difficult task, especially when N, K become large. 
Linally, computing A* as the fixed point of (3), (8) or (12) 
does not provide any insight into the optimal structure of A* 
and p*. To overcome these issues, we exploit the statistical 
distribution for h.ijk and and the large values of N, K to 
compute deterministic approximations of A* and p* [21]. Lor 

^This case is considered for example in [16] in which the SINR constraint 
is achieved by all UEs through an iterative procedure based on the bisection 
method. However, such a procedure can only be applied when the same SINR 
constraint is imposed to all UEs, which is not the case of the considered 
network. 


technical purposes, we assume the following grow rate of the 
system dimensions: 

Assumption 1. Ai V —)■ oo, 0 < liminf^v-j-oo AT/V < 
lim supjv_j.g^ K/N < oo. 

A known problem with the asymptotic analysis is that the 
target rates are not guaranteed to be achieved when N is 
finite and relatively small (e.g., [14]). This is because the 
approximation errors translate into fluctuations of the resulting 
SINR values. However, these errors vanish rapidly when N 
takes large yet finite values as it is envisioned for massive 
MIMO systems [2]. 

We further assume that only imperfect CSI is available at the 
BSs. Since the optimal linear precoder for the aforementioned 
network configurations is not known when only imperfect 
CSI is available, we overcome this issue by simply replacing 
the true channels with their estimates (which is an accurate 
procedure for good CSI quality). 


A. Coordinated Beamforming 

Let hijk be an estimate of hijk and assume, similar to [11] 
(among many others), that this can be modeled by the generic 
Gauss-Markov formulation: 


Mjk = ydljk 1 — T'fjk'^ljk +'Tijk(\ljt^ = ydijk^ljk (14) 


where c\ijk ~ CN{Q,1n) accounts for the channel estimation 
errors independent of wjjfc and 2,1 jk = — Ty-wijk + 

TijkCiijk- The parameter Tijk € [0,1] reflects the accuracy or 
quality of the channel estimate, i.e., Tijk = 0 for perfect CSI 
and Tijk — 1 for a channel estimate completely independent 
of the genuine channel. Replacing {h^Hj by {hja} into (3) 
yields 


\* _ 

^jk - 


(1 -I- l/^jk) 


“iife 


EEA}AHh|f,+iVI^. 

;=i i=i 


-1 


yj,k. (15) 


^jjk 


As mentioned earlier, the implicit formulation for A*^, prevents 
any insightful analysis of the system performance. By a large 
dimensional analysis, exploiting recent tools from random 
matrix theory (see notably [17]), we shall subsequently show 
that A*^ gets asymptotically close to an explicit deterministic 
quantity as N and K grow large as for Assumption 1, and 
that this quantity provides clear insight on the behavior of the 
precoder and the system as a whole. 

Lor technical reasons, the following reasonable assumption 
is imposed on the system settings. 


Assumption 2. The {dijk} and {ijk} satisfy 
limsupmax{(iyfe} < 00 and limsupmax{ 7 jfc} < 00 . 


A main technical result then lies in the following theorem: 


Theorem 1. Let Assumptions 1 and 2 hold. Then, 
maxjfe |A*^ — Aj^ ^1 0 almost surely with 


A 


{CoBF) _ 1 7jfc 


jk 


rjj djjk 


(16) 
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where {r]j} is the unique positive solution to the following set 
of equations 

-1 


Vj = 


L K dju 1 

dm m 

2—/ 2—d 1 I ^ 

1=1 i=i ^ + 7/* dn, m 


or, equivalently. 


L K dju Vj 


N ^ ^ 1 4- . iiiii ILi 

i=i i=i ^ + 7^^ dn, m 


Vj 


Vj- 


(17) 


(18) 


for A 


(CoBF) . 


•jk 


X 


•jk 


7A 


jk‘ 


L K 

m=^j--r Y. Y' 


1 

TV 


<iiH »)( 


1=1,l^j i=l 1 + 7H 


d-iii V] 
dm m 


where is defined as 


K 




l]k 


k=l 


-Ijk 


form as stated in the following corollary (similar results were 
obtained in [16] for a two-cell only network): 

Corollary 1. If a completely symmetric scenario is considered, 
then Vj r]j = 7} with 


L K 

1 1 
r? = 1- 

' AT 


'® dm 


AT 1 

1=1 i=l ^ 


'® dm 


( 21 ) 


Proof: The main difficulty lies in the implicit definition 
of the A*^’s in (15). A first step consists in heuristically 
discarding the implicit structure to retrieve the expression 


in explicit form. To proceed with an accurate 
proof, in Appendix A we follow similar steps as in [17] (in a 
completely different context though), by controlling the ratio 

■^(CoBF) 


Some important insights can be readily extracted from 
Theorem 1. To begin with, observe that the computation of 
ill (16) for cell j requires only the knowledge 
of the SINR constraints {juiVlji} and the average channel 
attenuations {djuiVlji} from all UEs to BS j. Differently 
from the fast fading channels {hju} required in (15), the 
latter can be accurately estimated as they change slowly with 
time (relative to the small-scale fading). Also, the Lagrange 
multiplier Xjk is known to act as a user priority parameter 
that implicitly determines how much interference a specific 
UE k in cell j may induce to the other UEs [7]. Erom (16), it 
turns out that is proportional to jjk and inversely 

proportional to djjk such that higher priority is given to 
UEs that require high performance or have weak propagation 
conditions. Moreover, observe that (18) can be equivalently 
rewritten as: 


Proof: The proof follows directly from Theorem 1 using 
djii — diji and yn — •yi Vj, 1. I 

An explicit form for {rjj} can also be obtained in the high 
SINR regime as stated below: 

Corollary 2. If yu grows large, then \/j 

lim\/id-yii=foo Vj = V with rj = 1 — KL/N provided 
that i -KL/N > 0. 

Proof: The proof follows directly from (17) or (18) of 
Theorem 1 observing that for given sets of {dju} and {dm} 

<i,-H Vj 


the ratio 


7H- 


' <illi VI 


1 as VT, i yii —>■ oo. 


We now proceed to computing the asymptotic powers 
{pfk°^^^} satisfying the SINR constraints in the large-(W, K) 
regime when imperfect CSI is available. To this end, we 
first compute the asymptotic values of the SINRs under the 
assumption that the transmit powers {pjk} are held fixed at a 

constant value and that gjk is replaced with gjk = 
such that 


SINR,fe = 


Pjk l^jjk^jkl 

N Hr,* 17 


K 

pj_ 

N 


L K 

+ E 

i=i,iy^j i=i 


where v*^ = 


(EtiEti 


K 


hjuh}}, +1 


= 1 N 


N 


Vj (19) 


( 20 ) 


then have the following result: 


-1 , 


^jjk- 


( 22 ) 


We 


Lemma 1. Under Assumptions 1 and 2, max^fc |SINRjfc— 

-(CoBF ), 


SINR 


jk 


► 0 almost surely with 


=Pjudjjk {I-rlk) 


L K 

i-^EE 


[Hi 


Hi ) 

"KuvNf 


l=lr=l (i+7h^^) 


Erom the above equation, it follows that r\j in cell j is such 
that r\j < and also it depends not only on its own cell 
requirements {yjk} through but also on all the other cells 
through the ratios {lii'^yf'^} and might be thought of as a 
cell priority parameter: higher priority is given to cell j if pj 
is small. In particular, djH /dm describes the relative strength 
of the interference received at UE i in cell I from BS j; for 
a given set of {yn} it is almost one for cell edge UEs of 
neighboring cells, while it is almost zero when cell I is very 
distant from BS j. In other words, higher priority is given 
to those cells that create high interference, as it should. To 
get further insights on this, please refer to the simple two-cell 
two-user network case study considered in Section III-E. 

If a completely symmetric scenario is considered, i.e., 
djii/dm — ^iji/^jji and yn — y^ \/j,l, the coefficients 
are all equal to a given p, which can be computed in explicit 


jiCoBF) 

^jk 


(23) 


where fkijk HXiPu) 


1 - 1-2 
-L Ljk 


Pljk = dijk 


(i + wfeS)' 


(24) 


Proof: Substituting X^^^ for the result follows 
as shown in Appendix B. ■ 

Eor notational convenience, let us now denote by b = 
[5i,..., 6^]^ the vector with entries 


1 ^ 

h = -y 

d AT Z^ 


7j7 




( 25 ) 
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The main result of this section unfolds from the previous 
lemma and provides the transmit power dedicated to each 
user and the minimal total transmit power to meet the SINR 
constraints. 


Theorem 2. Let F € be diagonal with entries 


L K 








and F € such that 


K 


A 1 IjkPljk 






(26) 


(27) 


where fdijk is defined in (24). If and only if 
limsupjf ||r~^F|| < 1, then under Assumptions 1 and 
2, the powers {pjk} required to meet the SINR constraints in 
the asymptotic regime are obtained as: 


(CoBF) 

Pjk 


Ijk 


Y. PljkPl 


-b(CoBF) 


di 


■jjk 


- 4 fc) 


L K 

i-^EE 


[Hi 


‘Lja2±' 

’ll . 


(28) 


Z—/ Z—/ / d-iu "ni \ ^ 

1=1 


collects the total transmit power of each BS with b defined as 
in (25). Moreover, the asymptotic total transmit power is given 
by 


piCoBF) ^ ^ji—(CoBF) 


= (F-F)’ 


b. 


(29) 


Eti I — [r']j j the strict inequality holding for 

at least one j. Therefore, it follows that if Vj € {1,..., L} 




K 


Tj k l^lj k 


< 


K 








(30) 


with the strict inequality holding for at least one j, then 
limsup;^ ||r“^F|| < 1. Assume now that the quality of the 
channel estimates is the same for all UEs, i.e., Tijk = tj VI, k, 
then from the above discussion and using (24) (after simple 
calculus) the maximum level of imperfect CSI in cell j is 
found to be; 


Lemma 2. Assume Tijk = Tj VZ, k. Then, for any given 
set of {yjk^ydjjk} and {dijk\ a feasible asymptotic power 
allocation exists in the CoBF case if Tj < Tj^^^'^ with 


(CoBF) 

EiHiax 


1 + ' 


K L 

Ijk 
‘jjk 


- 1/2 






=1 


L K 

i-^EE 


, 7H- 


•ijli Vi 


+7ii 


(31) 


1=1 i=l 


J 


To confirm the asymptotic optimality of the power alloca¬ 
tion provided by Theorem 2, we now provide the following 
corollary that proves that the uplink-downlink duality holds 
true also in the asymptotic regime when perfect knowledge of 
the channel is available: 


Proof: The proof is given in Appendix C and basi¬ 
cally proceeds as follows. The transmit powers are set to 
ensure that the SINR constraints are reached exactly in the 
asymptotic regime, that is such that = yjk 

(with defined in (23)). It then suffices to solve 

the implicit equation = yjk in the unknowns 

3"his equation turns out to unwrap as an explicit 
equation for the which are then readily obtained 

as in the theorem statement. The asymptotic approximation of 
the total transmit power easily follows taking into account that 
the transmit power of BS j is given by Yk=i ~ 

piCoBF) 

j " 

A close inspection of Theorem 2 reveals that the compu¬ 
tation of only requires knowledge of the system 

parameters {yu} and {rm} and the large scale channel fading 
components {dijk}- The latter change slowly in time com¬ 
pared to small-scale fading components and can be accurately 
estimated and possibly exchanged among coupled BSs with a 
reasonable effort. This is in sharp contrast to the finite system 
regime wherein the evaluation of through (4) requires 

some channel exchange procedure within the network at the 
rate of the fast fading channel evolution. 

From [22, Corollary 3], we have that limsup;^ ||r-iF|| < 
1 if none of the column sums of F ^F exceed unity and at 
least one is less than unity. This amounts to checking that 
Yi=i j / — 1 equivalently, (since F is diagonal) 


Corollary 3. Under Assumptions 1, 2 and 4, if perfect 
knowledge of the channel is available, i.e., Tijk = 0 yi,j,k, 
then the duality gap between (2) and (5) is zero. 

Proof: The proof is sketched in Appendix D. ■ 


B. Coordinated Multipoint Processing 

With a slight abuse of notation, let h.jk € be the 

estimate of the channel from BS j to user k given by 
hjfc = -f T kCjjk) = \/djkZjk where 

qjk ~ CA/'(0, Iat) ai^ Zjk - '^I'^jk + r^q^fc. Then, 

we may write = [h^^,..., ,kV as 

hfc = 0/^ -f Tfcqfc^ = ©/^z/c (32) 


with 0fc € £NLxnl defined as ©^ = 

diag{difc,...,dLfc} O Ia and q^ = , qE]^ and 

Wfc = [w^^,..., Replacing with hk in (32) yields 


At = 


(1 -f l/7fc) 


-1 


fKL 


(33) 


hf (^E Ath,hf+iVLIwij hk 


Vk G {1,2,..., KL}. Similar to the previous network config¬ 
uration, we shall require here the following technical setting. 


Assumption 3. The {djk} and {jk} satisfy 
limsupmaxj^fc djk < 00 and limsupmaxfc 7 fc < 00 . 
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Our first result in this setting is as follows: 

Theorem 3. Under Assumptions 1 and 3, max^ I At — 

-iCoMP ), 


0 almost surely with 

—(CoMP) _ 

= - 


Ik 


X Y. diktii 


(34) 


If and only /f limsup^^ ||Z|| < 1, then under Assumptions 1 
and 3, the powers {pk} required to meet the SINK constraints 
in the asymptotic regime are obtained as: 


(CoMP) _ Ik Cfc 

Pk 1 - 

^ ‘k 


flk 




1 - (1 +7fc)' 


(i+7fcr 


+ ^2 (41) 


where {pi\ is the unique positive solution to the following set 
of equations: 

_ / 1 dii 7i 


h-i = 


\NL^, L 

V .=1 ^ ^ 


7 = 1 


1 


■It 


+ 1 


V/. (35) 


1 ^ 

Efe = ^ ^ dikPi 
1=1 


NL ||9J||2 


KL 

Y WL 


P. |hfv*P 


where = 


(y- 


KL A* 


= 1 ■^^NL 


-1 , 


further, we call ej, = 

as e'f, = {Ikl ~ J)~^ Cfc where c'^, e is such that 


) ^KLk, 


''k\i 


Ik 


i Y dikPi 



where fl = computed as ft = 

ct^(I_r'l-Z)~^z and z G with Zk = t^Yp=i rr^^- 
Moreover, the asymptotic total transmit power is 

KL 


-^(CoMP) 
r 0-1 


= —y 

NT, 


NL 


(CoMP) 


(42) 


i=l 


Proof: The proof is given in Appendix E and proceeds 
as that of Theorem 1. ■ 

Unlike (16), in the CoMP configuration the Lagrange mul¬ 
tiplier of UE k is found to be inversely proportional to a 
weighted priority parameter given by 


Proof: The proof follows from Lemma 3 using the same 
arguments for proving Theorem 2. ■ 

As done before for both CoBE, we observe that if 


1 

fVL 


KL 


li 


T.-- 


^ik 


l-Tf 


1 - (1+7*)" 


(36) 




(1+7*)^ 


+ < 1 


(43) 


which basically takes into account the effort of each cell for 
jointly serving user k - see Section III-E for an illustrative 
example. Replacing with g^ = yjffp ||y*|p the SINR of 
user k takes the form 

SINRfe = 


and less than unity for at least one k, then limsup^^ ||Z|| < 1. 
Therefore, we have that: 

Lemma 4. For any given set of { 7 ^} and {dn}, a feasible 
power allocation exists for CoMP if Ti < with 


( 


(37) 


(CoMP) 


1 + 


7* 


z,max 
KL 

NL ^ ^ik 


V 


ef(l + 7*)%_^-,. / 

t NLIL 2^ ^ik 
‘ fc=l 


- 1/2 


(44) 


hfc. To proceed 
^ the vector obtained 


2 I T'y^^^iidikpi\ (38) 


Corollary 4. Under Assumptions 1 and 3, if perfect knowledge 
of the channel is available, i.e., = 0 Vfc, then the duality 

gap between (7) and (10) is zero. 

Proof: Despite being much more involved, the proof 
basically unfolds from the same arguments used for proving 
Corollary 2. ■ 


and J G has entries given by [J]j 

Our main results are then as follows: 

Lemma 3. Under Assumptions 1 and 3, max^ jSINRfc — 
SINR^^ ^ I —)■ 0 almost surely with 


SINR^ 


(CoMP) 


= Pk- 


l-'ffc 


I) 


(39) 


1 


li 


ik 


NLl-rf e? 


l-N 


1- (1 + 7*)" 


(1+7*)" 


(40) 


C. Single Cell Beamforming 

Replacing hijk with hijk into (12) yields: 


\ * _ 

Nk - 


(1 + l/7jfc) 


-1 


\i—l 


yj,k. (45) 


where — t ^fc[i (i+ 7 p)^] / 1 ^^kl ^ 

wnere ik — (l+fkf ( atl Z^*=lP*Tf 7 

e' = [ci, ■ • •, = (Ikl—J)~^ c where c G has 

elements [c] - = + Yf=i dupf. 

Proof: See Appendix E. ■ 

Theorem 4. Let Z G be such that 


Assumption 4. The {'yjk} satisfy limsupjy ^ Yf=i 1+2/■■ + 

When the Lagrange multipliers are computed as above, we 
have that: 

Theorem 5. Let Assumptions 1, 2 and 4 hold. Then, 
maxjfe |A*^ — Ajf ^ 0 almost surely with 


A 


(ScBF) _ 1 7jfc 


‘jk 


+ ^jjk 


(46) 


where 1 /, is given by ( 20 ). 



Proof: The proof is an easy extension of [12, Theorem 
1 ] derived for a single cell network. ■ 

As expected, the computation of for cell j 

requires only knowledge of the SINR constraints and average 
channel attenuations of its own UEs, i.e., {jjk;yk} and 
{djjk'yk}. Replacing in (45) with in (46), we 

then have the following result. 


Lemma 5. Under Assumptions 1, 2 and 4, max^^ jSINRjfc — 
I —)■ 0 almost surely with 


SINR 


{ScBF)_ 
■jk ~ 


K 


1—y 

1 N ^ 




Pjkdjjk (l - T^jk) 
jiScBF) ^ 

(47) 


where = Ylf=iC^ijk (^wYlf=iPii^ accounts for the 

interference with 

^ ~ ‘ (48) 

dijk j f I- 

Proof: The proof is the same as for Lemma 1, the main 
difference lying in the fact that V;* in (22) is now independent 
from h.ijk as shown in Section II-C. ■ 

Theorem 6. Let A G be a diagonal matrix with 

elements 



K 




''ft 


^ ^ (1 + lokY 

and U G such that its (j, l)th element is 

K 




Ijkdljk 


N ^ d,,k{l - Tjjk) 


(49) 


(50) 


where aijk is defined in (48). If and only if 
limsup;^ ||A~^U|| < 1, then under Assumptions 1, 2 
and 4, the powers {pjk} required to meet the SINR 
constraints in the asymptotic regime are: 


(ScBF) 

Pjk 



Z=1 


K 


1 _ 'y' 'ji 

N (1+7.0" 




(51) 


where = [pf 


fScBF)yj. 


py-—-Y=aYA-V)-^h 

collects the total transmit power of each BS. Moreover, the 
asymptotic total transmit power is 


P. 


iScBF) ^ ^yp(ScB+) _ ^2 iT 


= a^l^'(A-U)-'b. 


(52) 


Proof: The proof is similar to that of Theorem 2. ■ 

As seen, the computation of and 

quires the BSs to exchange almost the same system parameters 
of CoBF. Then, one might argue that ScBF has no potential 
advantage with respect to CoBF. The advantage comes from 
the observation that, unlike CoBF, the implementation of ScBF 
does not require knowledge of {hj^; V/ f j, i} at each cell j. 
Although it is true that this information could be potentially 


acquired at cell j, operating with a time-division duplex 
protocol, simply by reception of the pilot signals transmitted 
by all UEs, from a practical standpoint, this task would require 
a proper allocation of pilot sequences to avoid the so-called 
pilot contamination effect. In short, compared to CoBF, ScBF 
allows to simplify the channel estimation task. 

As done for CoBF and CoMP, we use [22, Corol¬ 
lary 3] to state that if ^ ^ < 1 - 


=1 d. 


1 

N 2^k^l 


Ink 


(1+7.. 


Vje{i,., 


L} or, equivalently. 


1 ^ 

f_ Y Tjfc 
^ fc=i ^ \ 




ijk I < 


(53) 


with the strict inequality holding for at least one j, then 
limsupj^ ||A^^U|| < 1. Imposing Tjjk = U Vfc and using 
(53), from the above condition it follows that (after simple 
calculus): 


set of {'yjk}, {djjk} and {d/jfc}, a feasible power allocation 

with 


Lemma 6. Assume that Tjjk = Tj Vfc. Then, for any given 
, {djjk} and {d, 
exists in the ScBF if Tj < rj 

(54) 


(ScBF) 




1/2 


V 

and Vj G {1,..., L}. 


K 

1 - - y" 

^ N 2^ 


Uk 


N f-' 1 + 7 .)= 

k—1 


J 


If we assume that jjk = yj Vfc, then from (53) the maximum 
value of 7 j in cell j as a function of the other system 
parameters is obtained as follows: 


Lemma 7. Assume 'jjk = 7j Vfc. Then, for any given set 
‘^f {pjk}> {djjk} and {dijk}, a feasible asymptotic power 
allocation exists for ScBF if jj < yj^max with 


41. 


Tj.niax — 


^ - 1 
N ^ 


2Aj 


1 + i /1 + 4 


A, 


A, 


dX - 1 

N ^, 


(55) 


and A^ = A J2k=i 


'=+'*),■)= y=Tiit3 ' 

-- 


1 — 7*, 


33 k 


The asymptotic optimality of the power allocation of The¬ 
orem 6 is confirmed as follows: 


Corollary 5. Under Assumptions 1, 2 and 4, if perfect 
knowledge of the channel is available, i.e., Tijk = 0 yi,j,k, 
then the duality gap between (11) and (13) is zero, i.e., 

j_ xr^K -(ScBF) _ j_ -TiScBF )_2 

N 2-^k—lPjk N 1 '^jk ^jk' 

Proof: The proof is sketched in Appendix G. ■ 


D. A simple case study 

We now consider a simple case study that allows to easily 
confirm some of the insights observed above (and to get 
further ones). In particular, we consider a two-cell network. 
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Fig. 1. Graphical representation of a simple two-cell network with one UE 
per cell such that dm = d 2 ii = d 22 i = d and di 2 i = oi-d with a < 1. 


i.e., L = 2, wherein a single UE is active in each cell K = 1 
with 7 ii = 722 = 7 - We assume also for simplicity that 
the network exhibits some symmetry (as shown in Fig. 1). 
In particular, we assume that UE in cell 1 is positioned at the 


same distance from BS 1 and 2 such that dm = d 2 ii = d 
whereas UE in cell 2 is such that 1^221 = d and di 2 i = ad with 
d being the resulting path loss and a < 1. With ScBF, we have 
that A'fwith c = <ri =^2 = 

In the CoBF case, it turns out that = —1 and 

’ Li r)i a 


A 


{CoBF) 


= with (from (18)) 


m 


m 


— ( ^ + 7 \ ^ 1 7 

^ Vl + 7 a+7/ +7 

1 [ 7 I - 1 7 

n[i + 1 1 + 7 / " ^ 1+7 


(56) 

(57) 


where we have defined rj = ? 72 /pi- Since a < 1, it easily 
follows (by contradiction) that rji must be larger than 772 such 
that All < A 21 . As intuitively expected, higher priority is given 
to cell 2 since it creates higher interference than cell 1 due to 
the shorter relative strength. As for the CoMP, it turns out that 
ei = d /2 (pi + 7 x 2 ) and 62 = d /2 + 7 x 2 ) from which one 

gets Ai = — 7 —^ and A 2 = —— /;■ Since a < 1, then 
X 2 > Ai. Therefore, higher priority is given to user 2, as it 
should since BS 1 is located further away from user 2. If a 
symmetric network is considered such that for example a = 1 , 
then it turns out that pi = 772 = P with 77 = 1 — and 

that ei = 62 = d /2 (/xi + 7 x 2 )- 

From Lemma 6 , the maximum value of imperfect CSI in 
the ScBF configuration is found to be; 


(ScBF) 
‘ l,max 


(ScBF) 
' 2,max 



(58) 


(59) 


from which, since a < 1 , it follows that (as it should be due 
to the higher interference experienced by cell 1 ) < 

■ Moreover, from Lemma 2 we obtain 


(CoBF) 

l,max 


1 + 


^7 


\ 


- 1/2 


1 -F 


^ , (?r+'r\ 

i+7 + -(tt^;/ 


(60) 


(CoBF) 
2,max 


- 1/2 


1 + 


:^7(l + a) 


_ J. ( _!_ 

N I 1+7 


1 ( 7 ^) +i 


(61) 


with ^ as it follows from standard analysis 

taking into account that 77 < 1 and a < 1. Finally, the values 
of and can be obtained from Lemma 4. 

The lack of explicit expressions for {ei} and {e^ /j} does not 
allow an easy comparison between the two values. However, 
numerical results can be used to confirm the intuition that 

(CoMP) (CoMP) 

U,max ^ ^2,max 


E. On the limiting case N ^ 00 and K/N ^ 0 

We now look at the limiting case in which N ^ 00 and 
K/N ^ 0. The following results are easily obtained from the 
asymptotic analysis above; 

Corollary 6. If N ^ 00 with K/N 0, then v*^ = hjk 
and 

J2 


(CoBF) _ 

Pjk - 
Also, we have that 


a 


Ijk 


NP. 


(CoBF) 


1 

-+ 0 such that 
2 


L K 


a 


= 2.2.— 


Hi 


Ui 


dm 


(62) 


(63) 


i=l i=l 

Proof: The proof follows from Theorems 1 and 2 ob¬ 
serving that if A^ —)■ 00 such that K/N —>■ 0, then rjj —)■ 1, 

^ 0 . 


-riCoBF) _ . a,x 

^jk — Ijk/djjk, 2.7=1 Z.i=l 


N 




Also, [r]jj- —^ 1 and 0 Vj, /. 

Corollary 7. If N ^ 00 with K/N —0, then and 


(CoMP) 

Pk 


7fc 


1 — L 

^iEdik 


(64) 


7=1 


, , T^(CoMP) 

Also, we have that Prp 


NP. 


■(CoMP) 


0 with 


KL 

E O Ik 

1 _ ,-2 L 

'==1 '' E dik 

7=1 


(65) 


Proof: The proof follows from Theorems 3 and 4. If 
N ^ 00 and K/N —>• 0, then )ii —1, ej. —>• ^ E2 i 


(CoMP) _ _ j ^ 


= 7 fe/(rE 7 =i'^ifc) and E^i 


A, 

Also, Ufe —)■ 0 since Zk —> 0. 

Corollary 8. If N ^ 00 with K/N 
down to CoBF. 


NL 


h,hf ^ 0. 


0, then ScBF boils 


The above corollaries state that, if N grows unbounded, 
the precoder for power minimization reduces to the maximum 
ratio transmit (MRT) scheme for all network configurations. 
Moreover, it turns out that the power required by all schemes 
to meet the constraints {jii-,yi,i} (or, equivalently, { 77 £;Vfc} 
for CoMP) vanishes inversely proportional to 1/N. Moreover, 
the following result holds true; 
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Lemma 8. If L > 2, N ^ oo and K is kept fixed, then 


NP. 


■jrr(CoMP) 


< NP. 


■j-{CoBF) 


= NP. 


-jrjiScBF) 


( 66 ) 


Proof: Extending to CoMP the double index no¬ 
tation used for CoBF, (65) becomes NP^°^^'^ = 

v=rrr~ 5 — from which the result follows 

since Y.j=i dju > dm when L > 2. ■ 

The above result states that the total power consumption 
decreases faster for CoMP than for CoBF or ScBF, meaning 
that a fully-cooperative system provides potential advantages 
for power saving. Also, ScBF performs as CoBF when the 
number N of antennas grows very large. Consider for example 
a system in which ^m = 7 and Tm = rVlf and, thus, jk = J 
and Tk = T Vfc. Assume also that dju = d In these 

circumstances, as N grows we obtain 




1 — d 


NPfl ' =K- ^ 


—(CoMP) 
T 


I — d 


(67) 

( 68 ) 


from which it follows that a power saving of 10 log L dB 
is achieved with CoMP in this particular setting. This can 
be potentially used to reduce the number of antennas by a 
factor L. The above results are only apparently in contrast to 
those in [23] wherein it is shown that a CoBF system provides 
considerable performance improvement (under a wide range of 
utility functions) as compared to CoMP. Indeed, the results of 
[23] are obtained under the assumption that the same number 
of spatial degrees of freedom per UE is provided by each 
network conhguration. This would amount to assuming that 
each BS in the CoBF and ScBF settings is equipped with 
NL antennas instead of N. Under this assumption, our results 
corroborate those in [23]. However, this would require to 
increase the number of antennas per BS by a factor L for 
CoBF and ScBF. 



Fig. 2. Average transmit power in Watt vs. target rate when L = 4, K = S 
and N = 32. 



Fig. 3. Average transmit power in Watt vs. target rate when L = 16, K = S 
and N = 32. 


IV. Numerical results 

Monte-Carlo (MC) simulations are now used to validate the 
above asymptotic analysis for a network with hnite size. The 
results are obtained for 1000 different channel realizations and 
UE distributions. We consider a multi-cell network composed 
of L square cells distributed in a square region of side length 
500 m. Following [12], the path loss function dijk is obtained 
as dijk = 2L^[1 + \\xijk\f/x'^) where xijk S is the 
position of user k in cell j with respect to BS Z, k > 2 
is the path loss exponent, t > 0 is some cut-off parameter 
and Lx is a constant that regulates the attenuation at distance 
X. We assume that k = 3.5 and = —86.5 dB [12]. 
Similarly, we have that djk = 2L^(l + with 

Xjk being the position of UE k with respect to BS j. The 
transmission bandwidth is VU = 10 MHz and the total noise 
power Wa^ is —104 dBm. Unless otherwise specified, in the 
subsequent simulations we assume that the same data rate must 
be guaranteed to each UE. Moreover, we assume that K = 8, 
N = 32 and impose the same quality of channel estimate for 
each UE, i.e., Tm = t yi,i. 


Fig. 2 illustrates the average transmit power in Watt vs. 
target rate in bps/Hz/UE when L = 4. Markers are obtained 
using the asymptotic analysis whereas the error bars indicate 
the standard deviation of the MC results. Clearly, = 0 
corresponds to the perfect CSI case. As expected, the higher 
power consumption is required by ScBF due to the lack of 
cooperation among BSs. Compared to CoMP, a slight increase 
of power is required by CoBF. As seen, the approximation lies 
roughly within one standard deviation of the MC simulations 
and thus we may conclude that the large system analysis is 
accurate even for networks of finite size. 

We now investigate the performance of the different 
schemes when the network becomes denser. To this end. Fig. 
3 plots the average transmit power in Watt vs. target rate when 
L = 16 such that the total number of UEs in the network is 
KL = 128. Similar conclusions as for Fig. 2 can be drawn 
with the only difference that the average transmit power for 
target rates up to 2.5 bps/Hz/UE is smaller for all schemes 
due to the shorter distances of UEs from their serving BSs. A 
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Fig. 4. Average transmit power in Watt vs. when L = K = 8, N = 32 
and r = 2 [bps/Hz/UE]. 



Fig. 5. Average transmit power in Watt vs. when L = 16, K = 8, N = 
32 and r = 2 [bps/FIz/UE]. 


larger gap between CoBF and CoMP is observed compared 
to the results of Fig. 2. This means that full cooperation 
is benehcial as the network becomes denser. This is true 
especially when imperfect CSI is available and high target 
rates must be ensured as it can be deducted by the rapid 
increase of the transmit power for CoBF when =0.1 and 
r > 3. Moreover, it can be seen that ScBF cannot support 
target rates beyond 3.5 and 3 [bps/Hz/UE] for = 0 and 
= 0.1, respectively. 

Figs. 4 and 5 plot the average transmit power in Watt 
vs. accuracy of CSI when the target rate for all UEs is 2 
bps/Hz/UE. As expected, the average transmit power increases 
as T becomes larger and the CoMP technology provides more 
robustness to imperfect CSI compared to CoBE and ScBE 
This holds true also for larger values of target rates. 

Eig. 6 plots the average transmit power in Watt vs. N when 
K — 8 and L = 16. The label I-CSI refers to the imperfect 
CSI case. The values of ;V/,/} and {rii-,yi,i} are ran¬ 
domly taken in the intervals [0.01,0.1] and [1,3] bps/Hz/UE, 
respectively. As expected, the average transmit power of all 



Fig. 6. Average transmit power in Watt vs. N when L = 16, K = 8 and 
the rates are randomly taken in the interval [1, 3]. 


schemes decreases as N grows large but, for a given N, CoMP 
requires much less power compared to CoBE and ScBF. The 
latter perform almost the same as N is large enough. 

V. Conclusions 

In this work, we analyzed the structure of the optimal 
linear precoder for minimizing the total transmit power when 
different degrees of cooperation among BSs are considered; 
single cell processing, coordinated beamforming and coordi¬ 
nated multi-point processing. Stating and proving new results 
from large-scale random matrix theory allowed us to give con¬ 
cise approximations of the Lagrange multipliers, the powers 
needed to ensure target rates and the total transmit power. 
Such approximations turned out to depend only on the long¬ 
term channel attenuations of the UEs, the relative strength of 
interference among BSs, the target rates and the quality of 
the channel estimates. Numerical results indicated that these 
approximations are very accurate even for small system dimen¬ 
sions. Applied to practical networks, such results may lead to 
important insights into the system behavior, especially with 
respect to target rates, CSI quality and induced interference. 
Moreover, they can be used to simulate the network behavior 
without to carry out extensive Monte-Carlo simulations. It 
is worth observing that the asymptotic analysis provided in 
this article could potentially be extended to more advanced 
channel models to account for example for the correlation 
among BS antennas [11] and/or for the presence of line of 
sight components [24], [25]. All this is left for future work. 

Appendix A 
Proof of Theorem 1 

Let us first provide some intuition on the main result. To 
this end, apply the matrix inversion lemma [26] to rewrite (15) 
as follows 

¥ = { 4 E V.KjHk'!. + In') (W) 
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Assume erroneously for a moment that the scalars {A^j} are 
given and independent from the channel vectors {h^zi}. Then, 
using classical random matrix theory results such as those in 
[11, Lemmas 4 and 5] yields 


'jk \ 1=1 i=l 


-1 


Ijk 

X 


H 

■jli 


Using [11, Theorem 1] we have that 

L K 


t-N 


-1 


\ / = 1 2=1 


H 

jli 


i-N 




(70) 


(71) 


where the coefficients {rj*} are solutions of the following 
system of equations 


^7 = 


L K 

■EE 




-1 


1 + Xl^djuT]* 


(72) 


From the above discussion, we may then expect the terms A*^ 
to be all close to A*i, = -k- for N, K large enough. This 

Jfv rij djjk 

statement is made rigorous in the following. To this end, let 
us define 


Cjk — 


-^(CoBF) 

■^jk 


Ijk 1 
djjkVj A*fc 


(73) 


where the coefficients { 77 ^} are such that: 




-(CoBF) 

_ kljli 

N , -riCoBF) 

z=i i=i 1 + Xu djuilj 


/ . L K 

- vEE 


+ 1 


From (69), using (73) one gets 


(74) 


IjkCjk _ J_ H 
Xjkdjjk ~ N 


jy + Itv] Zjj/c (75) 


where the superscript is omitted for simplicity. As¬ 

sume that {cjk} are well defined, positive and such that 
0 < Cii < C12 < • • • < CiK < C21 < C22 < • • • < C2K < 
••• E cli < cl 2 < ••• < clk- Then, using monotonicity 
arguments, from (75) it follows that for j = L and k = K we 
have 

dLLK ~ 'N^kLK^k^icLK)zLLK- (76) 


with Bfe(cLK) — d- Clk'^N- 

Assume now that clk is infinitely often larger than 1 + f with 
f > 0 some positive value [17]. Let us restrict ourselves to 
such a subsequence. From (76), using monotonicity arguments 
we obtain 


Ilk 

XhKdLLK 




(77) 


Applying [11, Theorem 1] one gets 

■^^LLK^k^i^ + £)zllk ^ eL{e) (78) 


with Cj {£) being the unique positive solution to 


L K 




Xlidjli 




+ 1 +. 


(79) 


From (78), recalling (77) yields limi<-_>oo inf eL(f) > 
T^dXLLK' Using the fact that ez,(0) = rju and that euiH) is 
a decreasing function of £, it can be proved [17] that for any 
f > 0 limif_>.oo supeL(f) < This however goes 

against the former condition and creates a contradiction on the 
initial hypothesis that clk > 1 + f infinitely often. Therefore, 
we must admit that clk < 1 + f for all large values of K. 
Reverting all inequalities and using similar arguments yields 
Cli > 1 — ^ for all large values of K. Putting all these results 
together yields 1 — f < cli < cl 2 E • • • < clk < 1 -f f from 
which we may write maxfc=i 2 ,...,ri: \cLk — 1 | < ^ for all large 
values of K [17]. Taking a countable sequence of £ going to 
zero, we eventually obtain max^^i 2 ,...,i<' |cz,fc — 1 | —)■ 0 from 
which using (73) and assuming that lim;f_).oo sup 7^'° < oo 

—{CoBF) 

it follows that maxfe=i_ 2 ,...,rf \7^Lk ~ ^Lk \ 0 "'hh 

= 'jz^LLk' Tollowing the same steps for j = 
1 , 2 ,..., L — 1 completes the proof. 


Appendix B 
Proof of Lemma 1 

To begin with, rewrite the numerator of the SINR in (22) 
as follows 


Pjk 




N 




2 “ Pjk 




Afh 

N jjk - 


(80) 


-J-jjk 

where we 


with Aj = 

have replaced = Xji for notational simplicity. Ap¬ 

plying the matrix inversion lemma first [26] and then using 
(14) we have that the numerator in the right-hand-side of the 
above equation takes the form: 


—1 

N 


^jjk-^j^jjk — 


1 _ t-2 A 

^ 'jjk’^jjk-^j 


[jk] 1 


^jjk 


l + ^A,zchj5,Aflh 


jjk 


1 + '£ 


where Aj = 


( N Ei=l ^li^jli^fli + In 


jjk 

-1 


(81) 


and 


: [ok] _ ( 1 


= ^ Az^hjz^hjf,-f Izv 


(l,i)^{j,k) 


(82) 


According to [11, Lemmas 4 and 5] and the results 
of Appendix A as well as Theorem 1 it follows that 

vjjk X djjkijj, ^Xjkii^jL-^^j ^^jjk ^ Ij and 


2-h^ 

A 

N’^jjk^j 


\jk] 


^jjk 
1 

'pf^jjk 


A,h 


0. Therefore, we have that 

djjkLjj 


j’^jjk 


1 — t; 


jjk 2^ _|_ 


Ijk 


(83) 
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We proceed to computing the deterministic approximation of 
■^h.^jf.A'jh.jjk in the denominator of the right-hand-side of 
(80). Observe that applying the matrix inversion lemma [26] 
and using [11, Lemma 4] together with the results of Theorem 
1 (such that -fjk = ^jkdjjkVj) one obtain 

(l+7jfc)2 

with itr(A2) = A^tr(Aj^ - 2;lAr)-i|^=o [2, Theorem 4], 
Using similar arguments as those of Appendix A, we have that 

^tr (ajI - zIat) x?7j (z) (85) 

with ,y,(z) = U Eti Eti +1 - By 

differentiating along z, we get 


Using the resolvent identity Aj — Aj * = —Aj(Aj ^ — 
A^^^ )Ay^^ [26] and observing that 


A7^ - At'"'" 


jY (^'^ohjjkhjj/. + ciqjjkC[jj/. 


+C2hjjkqfjk + C2q]jk^fjk 


with Co = 1 - Cl = and C 2 = ^fc, from 

(89) one gets 

— A ( Ah- 

~ — AjfeCo —hjjfcAjhjjfc —— 


( L K E.-fLii'] 




33^^ J ]\J- '^33f^ 3 33'^ 

— ^jkC2-|^^jjk■^j^.j3k '^^jjk^j^jjk 


+ 1 ( 86 ) 


from which setting z = 0 and using simple calculus one obtain 
(omitting the functional dependence from z = 0) 


where we have defined Bj = A^:^ ' ' )Aj. 

Using [11, Lemma 7] we obtain that x 

u’{l + \jkCiu) J_iH A Vi ^ u{l+\jkCiu) J_iH An ^ 

1+A,.„ ’ AT^iife^ULjfc - l+\,.u ’ N’^^Jjk^J^Jjk - 


L K A,.iil±!LLY 

-r_ L V ' <iiii ) 

N Ai 2^ ( dju vj 3 ^ 
m=ii=i Y+'riiiqi--zi) 


_l-\-XjkU ^ 33^^ 3 33 l-\-XjkU 

— A,-fcC 2 M^ . H R Vl , ^ —\jk.C 2 'U,u' 

1+YkU N^nk^J^Jjk - 

defined, for notational convenience. 


where we have 




Therefore, itr(A^) x ry' such that j^h.f-f,A^h.jjk x 
<L3fc>i3 . Putting this result together with (83) we have that: 

(-l-+7ifc j 

Pjk ^ 

N ||v*J|2 

, 2 J, 1T-T' 'l 

Pjkdjjk {i Tjjk) 1 ?yz2z2/ , \2 ' 

We now deal with the intracell interference term in the 
denominator of (22), which can be rewritten as 


lHf]pf]Hf]"jA,h,,fc (88) 


with = [hjji,...,hjj(k-i),'^jj{k+i),---hjjK] G 


u' = ^tr + o(l). (93) 

Using the above results, (91) can be approximated as 

4i>s-.'A ( 4®““'pf'fif'") A,h,j* * 


^ Vi^ A f ^ 


m'( 1 -f XjkT^jkU) 


1 - 

~33k 


1 + ^jkU ^ 

or, equivalently, after simple calculus 


■ik , - ,2'’ 

(1 + ^jkU) 


K 

E 

P3^ 

N lie* Ip 




^ Vi^ A ■ 

^iijjkAj 


—h^ A - {1 A h 1. 

JY ’^33 fc TV J ^J^nk 

^-'^33k\^- i^+^jkUy 

- - -=-2- - 

(l + Xjkuj 


c— and Pfl ^ diag{^,...,^^£^, 

1 112 ) • • •) 1 ifif 112 }• In order to eliminate the depen- 

N I |■^J(fc+l) II iV I I I 


dence between hjjk and Aj, rewrite (88) as 
—h'^, A (A h i. 


_ J_Vi^ A^f^l (}_-\J(\jk]-p\jk]-^[3kY 


Ajhjjk + 


Notice that u in (92) is such that u x -^j^rjj (as it follows 
from Theorem 1) while the deterministic approximation of u' 
is computed as follows. Observe that 

N J f ^ J ) 

1 V n-.MhlYJl _ ^ (96) 

N ^ Av,H a2L.. ^ 

2=1 ,N 

from which using [11, Lemmas 4, 5 and 6] one gets 


y^k (a, - Af 1)( Inf Ipf iHf 1") A,h,,,. (89) 


n'-'-331 3 Uj* 


(V A2h -1. 
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Plugging the above result into (93) produces u' 


K 


djjk{-^ ^ Pji). Therefore, recalling (46) we have that 




—hY.uAi ( 


N 


iife J 




N ^ 


A,h 




Ijjk - 




K 


(1 + l]k) 


N 


Y.pi^ 


(98) 


i=l 


We are left with computing the deterministic equivalent of the 
intercell interference generated to user k in cell j by all other 
UEs in cell I ^ j: 


K 

E 


PH 1 


= -Kk^i ( -H,P,Hf ) Aihijk (99) 


N' 


H 


N 


with H/ = 
diag{^ 


IIK\ 


^NxK 


Pii 


PlK 


and P; = 
r}. Mimicking the same steps of 


above one gets 
1 




N' 


Hjk^l 


H 


N 




A,h 


1 - T, 


Ijk 


1 — (l + Xjku) 


(l + Xjku) 


( 100 ) 


where we have that u = ^lj|^tr(Ap*^) + o(l) and u' = 
l^tr(-^P/HpAp^^ H/) + o(l). Observe now that u x 
dijkTji and u' x Putting these results together 

and recalling (16) we eventually obtain 




1 - 


1 

A 
1 - 


H,P,Hf) AM,k 


^Ijk ■ 




K 


/ 1 • ^101) 
Putting all the above results together completes the proof. 


Appendix C 
Proof of Theorem 2 

From (23), omitting the superscript iC:°BF) jj; follows that 
the power p^f. such that SlNRjfc = is obtained as; 


Pjk 


Ijk 


’i (l - 4 *) 


{Ijk + CT^) ( 102 ) 


with rf^ given by (87). Plugging the above result into Pj = 

W Y.k=i Pjk an^observing that Ijk = Ya^i PijkPi it follows 
that the values {Pj} must satisfy the following set of equations 

Vj: 


K 


pA = 1V 


Ijk 




^ djjkil - T^fc) 


'^Pljkpl 


= E^'[Fk 


1 + cr bj 


(103) 


1=1 


where bj and [F]j; are given by (25) and (27). Rewriting the 
above set of equations in matrix form yields (F — F) P = a'^h 
with F diagonal and given by (26). Using [22, Theorem 2.1], 
a necessary and sufficient condition for a solution P > 0 to 
exist is that the spectral radius of F^^F is smaller than 1. The 
asymptotic approximation of the total transmit power easily 
follows taking into account that the transmit power of BS j is 
given by ^Yk=iPjk = Pj- This completes the proof. 


Appendix D 
Prooe of Corollary 3 

We begin by rewriting (F — F) 

A”^ (F — F) = (j2a where A G is diagonal 

with [A\ij = rii and a G with ai = 

Then, we have that 


Fa-‘( r-F)= 1 mi) 

1=1 i=l 

Observe now that the jlh entry of the row vector 
1^A~^ (F — F) is equal to unity Vj: 


i"A-^(r-F) 

(a) Pj 


dj 

Vj 




N 


Pj 


1=1 i=l 


1 +7H^ — 

, dm rji 


1 L K dju rij 

1 -^EE ' “ dm rii 


1 = 1 i=l ^ ^ 7'* dm Vi J 


(.)l 


where (a) follows from (24) - (27). The expression in (b) is 
obtained taking into account (87) whereas (c) follows from 
(18). Plugging the above result into (104) yields the statement 
of the corollary since the transmit power of BS j is given by 

piCoBF) _ j_ (CoBF) 

j N Z-^k—l^jk 


Appendix E 
Prooe of Theorem 3 

As done for proving Theorem 1. Assume erroneously for 
a moment that {A*} are given and independent from {hi}. 
Then, applying [11, Theorem 1] to (33) we conjecture that 

j(CaMP) 

the ratio Vk = is equal to unity with tk 

being the solution of the following set of equations: 

1 ' 1 

efc = - ^ dik^i = (GfcT) (105) 


KL 


with T = (^ x; 




2=1 




+ ^nl) From (33), using the 


IkTk = —2,% 0 ^ 


matrix inversion lemma and (32) one gets 
NL 

1,2 




J-ATL 


-1 




(106) 


where the superscript is omitted for simplicity. As¬ 

sume that 0 < ri < r 2 < ■ ■ ■ < tkl and that tkl is infinitely 
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often larger than 1 + f. Let us restrict ourselves to such a 
subsequence. From (106), replacing with tkl and using 
monotonicity arguments we obtain 

IKL < (107) 

with 

+ (1 +^)ItVL. 

i^k 

Applying [11, Theorem 1] yields 

X Pkl{£) (108) 

with pkl{£) being the unique positive solution to pk{£) = 
j^tr(0feT(f)) with 


KL 


The proof proceeds as in Appendix A. Omitting the details, 
we eventually obtain that 1 — £ < tkl < 1 + £ from 
which we may write max|rrfL — 1| < for all large values 
of KL [17]. Taking a countable sequence of £ going to 
zero yields maxjrifi — 1| —> 0 from which using tkl = 
and assuming Iitukl^oo sup l < oo we obtain 

u\&yi\X*j^]^ — \kl\ —> 0 with \kl = Following the 

same steps for A: = l,...,ArL — 1 yields the desired result. 


Appendix F 
Proof of Lemma 3 

We begin by computing the deterministic equivalent of 
the interference term of (37). This will be instrumental to 
obtain the deterministic equivalent of the useful signal power. 
We start rewriting the denominator of SINR^ in (37) as 

= (TOEt^A.h,hf + HW ^ 

[hi,...,hfe_i,hfc+i,...,hK] e C^x(A-l) p[k] A 

Pi Pk-i Pk+i 

oiagtYTTiivJIP J J Ji2j j|2) ■ • ■ 7 /• 


E KL in^ Vir 

with A 


NL I 


N L I 


Then, mimicking the derivations for CoBF, we have that 


1 


hfA 


1 

Wl^ 


.gWpWHW" ] Ah 


1-^fe 


1 — (l + Afcu) 


(l + Xku) 

where u = ■^tr(0fcAl^l) + 

j^tr(j^pWH['=lAW0feA['=lHW") 


0 ( 1 ) 

+ 

t-i 


and 

0 ( 1 ) 


(109) 


with 


Observe that 


At'^l = (i^EtM^fcEh,hf + 
u ^ €k (as it follows from Theorem 3) whereas the 
deterministic equivalent of u' is obtained as follows. Rewrite 
j^tr(j^pWH['=lAW0feA['=lHW"' 


as 


—iri J-PWhWa['=10,,AWhW"^ = 


NL 


— E 

NT. 


\NL 

KL 


NL 


Pi 


NL 


hf A0frAh, 




— h^ A2h 

NL’^i ^ 


( 110 ) 


Then, observe that 
1 


NL 


h,^A0t.Ah, 


(q 


0,A0t^A 


1 + ^A,tr 


(®^A) 


2 • 


( 111 ) 


To compute ■j^tr(0iA0fcA), 
l^tr(0,A0fcA) = j^£tr(0,(A-i - z& 

Observe now that [11, Theorem 1] 


we may write 

T ) ) U—0 


^tr (©. (a-‘ - ze, 


Tik{z) 


(112) 


where eik{z) is given by eik{z) = -j^tr {&iTk{z)) and 
Tk{z) is computed as 


Tk{z) = 


KL 

■E 




-1 


AT 1 + A„e„fc(z) 


lNL-Z®k . (113) 


By differentiating along z, we have e'j,(z) = T^tr (0iT).(z)) 


is given by T).(z) = 


where T',(z) ^ 

Tfe(2)(l^EEi (i+xblloir ^ ^ ^ 

yields 


4(O) = ^tr(0,n(O)) 


where 

4(0) =T 


_ rj, / 44(0)®?! 


E 




+ ®fc T. 


(114) 


(115) 


In writing the above result, we have taken into account that 
T = Tfc(O), e_„fe(0) = j^tr(0„Tfc(O)) = 4:1^ (®"T) = 
e„ and 7 „ = A„e„. Plugging (115) into (114) and neglecting 
the functional dependence from z = 0, e). = ..., 

is found as the unique solution of: 


,^^^i(1+7E 


®r 


— (®iT0fcT) 

I ^ t.' 


E 


n^nk 


41 (1 + Lnf NL 


tr(0iT0„T) . 


(116) 


Observing that ;4f''(®i'F0fcT) = j^{Yzi=i(^iidikp‘f) we 
may rewrite the above system of equations in compact form 
as e'f. = {Ikl ~ J)”^Cfc c). and J being defined in 

the text (see (38)). On the basis of above results, using 
4^Aitr(0iA) X 7 i, we eventually obtain that 


1 


A 01. Ah, 


^ik 


NL ‘ ■" " (l + 7i) 

Following similar arguments of above yields 


2 ■ 


1 

El' 


h^A^h, 


(1+7?)^ 


(117) 


(118) 
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with e' = [e'l, ..., = {Ikl — J) ^ c where c G 

L 

has elements [c]j = Putting the results in (117) 


1=1 

and (118) together, we have that 


1 

Wl 


^ ^hfA0feAh, 


1 

Wl 


KL 


ik 




(119) 


The deterministic equivalent of the numerator of the SINR in 
(37) is now easily obtained as 


Pk 


iTOhfvfcp 


■■Pk{l-rlf-^ 


( 120 ) 


since ^hfvfc x and j^hfA^hfc x 

4 

(l+7fc)"' 


Appendix G 
Proof of Coroflary 5 

To begin with, we use (46) and (48) to rewrite (51) as 
follows 


Pjk — 


^jjkPj + ^jk 

K 7? 


1 - — y" 

N 1^1= 


( 121 ) 


(1+7)7? 


where = J2i=i dijkPi+cr"^ and the superscript (■5°^^) 
has been omitted for simplicity. Plugging the above expression 
into Pj = ^Wli^=iPjk and solving with respect to Pj one 
obtain 


P^ = 


'1 - — 

N 1^1= 


7,-i 


N ^i=l ( 1 + 7 ,-,) 
? 


“ TY X] WkOijjk j ^ X 


f2 

^jk- 


( 122 ) 


fc=i 


k=l 


Replacing ajjk with its expression in (48) and recalling (46), 
it follows that (after simple calculus) 


1 _ A 

N Z^i=l 


"y ji 


AT X — 1- 


K 




N 


(123) 


k=l 


Plugging the above result into (122) completes the proof. 
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